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مقاله علمی- پژوهشی 
مقایسه روش‌های 181,۲) و بیزی در برآورد ارزش‌های اصلاحی ژنومی با معماری‌های 
۰ حلم ژنتیک 


رضا سید شریفی 3 فاطمه علاء نوشهر نعمت هدایت ایوریق ‏ حمال سیف دواتی! 


تاریخ دریافت: ۱۳۹۸/۰۳/۰۱ 


چکیده 


اتتخاب وتومی با | ستفاده ار 7 شانگرهای با تراک بالاه به ضومن 9۱1۳هایی که کل زنوم را بو هش می‌ذهند و اغلب درعدم صاذل پیو. شتگی با 
011های مجاور خود قرار دارنده ارزش ژنتیکی کل را پیش‌بینی می‌کند. این پژوهش با هدف بررسی تأثیر عواملی چون تراکم نشانگرهاء تعداد 011 
ورائت‌پذیری صفت و نوع توزیع اثر 0011 بر صحت برآورد ارزش‌های | صلاحی ژنومی با ا ستفاده از داده‌های شبیه سازی در گو سفند انحام گرفت. به 
همین منظور ژنومی متشکل از سه کروموزوم. هر یک به طول ۱۰۰ سانتی‌مورگان با سه مقدار ورائت‌پذیری برای صفت مورد بررسی (۰/۱» ۰/۲ و ۰/۵) و 
سه پنل ذشانگری (۸۵۰۰, ۱۰۰۰ و ۱۵۰۰) در سه سطح تعداد 011 (۵۰ ۱۰۰ و ۱۵۰) با دو اثر توزیع یکنواخت و گاما برای م0011 شبیه سازی شدند. 
صحت ارزش‌های | صالاحی ژنومی برآورد شده با | ستفاده از پنج روش نآ بیزش بیزظ بیزن) و بیز 1650 مورد مقایسه قرار گرفتند. نتایج 
این تحقیق نشان داد که هر چه تراکم نشانگر و ورافت‌پذیری صفت افزایش يافته و تعداد 4071 موثر بر صفت کمتر باشده صحت ارزش اصلاحی برآورد 
شده بالاتر خواهد بود. در بین روش‌های آماری زمانی که تعداد ,011 موّثر بر صفت پایین است و توزیع اثر 0011 گاما در نظر گرفته شد» پیش‌بینی 


ارزش‌های اصلاحی ژنومی با روش بیز1 عملکرد بهتری داشت. 


واژه‌های کلیدی: لانانان) بيزش بیزظ بیزن) بیز۸8۹0] ارزش اصللاحی ژنومی. 


و 


معقدمه 


در روش‌های کلاسیک اصلاح دام انتخاب برای صفات مهم 
اقتصادی با استفاده از اطلاعات شجره‌ای به همراه رکوردهای فنوتیپی 
خود فر آیندگان (فرزندان و نوه‌ها) و گذشتگان (والدین و سایر افراد 
خویشاوند نسل‌های قبلانجام گرفته و بهترین پیش‌بینی ناریب خطی 
از ارزش‌های اصلاحی (*01.[1" فرد حاصل می‌شود ٩(‏ با پیشرفت 
ژنتیک مولکولی استفاده از اطلاعات نشانگرهای ۸ در برنامه‌های 
اصلاح‌نژاده تحت عنوان انتخاب به کمک نشانگر (۷۸۹)" پیشنهاد شد 
(۳). بعد از آن شالوده انتخاب بر اساس کل ژنوم در مقاله‌ای که توسط 
موویسن و همکاران (۱۰) انتشار یافت» ارائه شد. در این روش می‌توان 
ارزش‌های اصلاحی ژنومی تمام افراد را با استفاده از یک مدل خطی 
ب‌صورت تابعیت فنوتیپ‌ها از تمام نشانگرهای متراکم که کل ژنوم را 


۱-دانشیار گروه علوم دامی دانشگاه محقق اردبیلی» اردییل» ایران. 

۲- دانش آموخته دکتری ژنتیک و اصلاح دام دانشگاه تبریز» تبریزه ایران. 

(#-ایمیل نویسنده مسئول: حمه.0مطه60) 1 تمطع0ع 56 _۲62۵) 
12577121262 1۱۳۵:۱022067 


پوشش می‌دهند با صحت بالا پیش‌بینی کرد. مهم‌ترین مرحله این روش 
برآورد اثرات آللی می‌باشد. زیرا هر چه اثرات .001 * به نحو صحیح‌تری 
برآورد شون ارزش‌های اصلاحی پیش‌بینی شده بیشتر معرف ظرفیت 
ژنتیکی حیوان بوده و رتبه بندی حیوانات به‌طور صحیح‌تری انجام خواهد 
شد و در نتیجه پیشرفت ژنتیکی و پاسخ به انتخاب بیشتری مورد انتظار 
خواهد بود. 

عوامل مختلفی می‌توانند صحت ارزش‌های اصلاحی ژنومی را 
تحت تأثیر قرار دهند. این عوامل شسامل مدل آماری مورد استفاده, 
تراکم نشسانگرهاء ورات پذیری صفت. تعداد داده‌های فنوتیبی در 
جمعیت مرجع (گروهی که دارای اطلاعات فنوتیپی بوده و از اطلاعات 
آنها جهت برآورد اثرات نشانگری استفاده می‌شود)» فاصله زمانی (تعداد 
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۳۳۲ نشریه پژوهشهای علوم دامی ایران جلد ۰۱۲ شماره ۲ تابستان ۱۳۹۹ 


نسل) بین جمعیت مرجع و جمعیت تأیید (گروهی که فاقد اطلاعات 
فنوتیبی بوده و ارزش اصلاحی ژنومی آنه تها بر اساس اطلاعات 
ژنوتیپی برآورد می‌شود) می‌باشند (۱۸). 

تحقیقات ذ شان داده ا ست که با افزايش تراکم ذ شانگر در سطح 
ژنوم» برآوردهایی با صحت بالاتر حا صل خواهد شد. | صولا افزایش 
تعداد نشانگرها در سطح ژنوم احتمال پیوستگی بین نشانگر و ب01 را 
افزايش می‌دهد که این مسئله باعث افزايش صحت برآوردها می شود 
(۲). همچنین در ورائت‌پذیری بالات ارزش‌های اصلاحی افراد سبهم 
بیش‌تری در بروز رکوردهای فنوتیپ داشسته و لذا در برآورد اثرات 
نشانگرها با استفاده از اطلاعات فنوتیبی صحت بیشتری مشاهده 
خواهد شد (۱۰). 

امروزه تلاش پژوهشگران بر آن است تا با شناسایی عوامل موثر 
بر صحت براوردهای ارزش‌های | صلاحی ژنومی بهترین بهره‌برداری 
از داده‌های موجود انجام شود. از آنجا که توالی‌یابی کل ژنوم و در 
اختیار دا شتن ارزش ژنتیکی افراد به خصوص در گوسفند به دلیل نوع 
استفاده از شبیه‌سازی روشی منطقی برای چنین موضوعی خواهد بود. 
ژنومی صفات پیوسته در معماری‌های مختلف ژنتیکی شامل توزیع‌های 
آ ماری مختلف اثرات ژنی» تعداد متفاوت ب011» سسطوح متفاوت 
«لنانان و بیزی شسامل بيزش بیزظ بیز ) و بیز 1۸850 انجام 
گر فت. تا علاوه بر م‌قاسه نحوه عملکرد روش های مختلف در 
معماری‌های مختلف ژنتیکی, بتوان میزان تراکم نشانگر و تعداد ,0011 
کارآمدی را برای برنامه‌های شبیه سازی جمعیت‌های گو سفند معرفی 


کرد. 


مواد و روش‌ها 
شبیه‌سازی جمعیت 

جهت ایجاد عدم تعادل پیوستگی مورد نظر بین نش‌انگرها و 
ب۵11ها از روش ویلامسن و همکاران (۱۸) استفاده شد. در این روش 
از عامل کو چک بودن جمءعیت و رانش, جهت ایجاد عدم تعادل 
پیوستگی استفاده می‌شود. به همین منظور جمعیت پایه کوچکی (00) 
با ۱۰۰ رس دام شامل ۵۰ رس نر و ۵۰ رأس ماده در نظر گرفته شد. 
فراوانی آللی اولیه برای چند شکلی‌های تک نوکلتوتیدی در سل پایه 
۵ (حاظ گردید. گامت‌های والدی بر اساس فرض عدم تعادل 
پیو ستگی و با ا ستفاده از روش تابع مکان‌ابی هالدان (۸) شبیه سازی 
شدند و سپس گامت‌های نوترکیب ایجاد شده به طور تصادفی انتخاب 
و آميخته شدند تا فرد جدیدی از نسل ) ایجاد شسود. روش مورد 
استفاده تا ایجاد نسل ۵۰ام پیگیری شد. بعد از ۵۰ نسل آمیزش 


تصادفی اندازه جمعیت به تعداد ۱۰۰۰ فرد (۵۰۰ نر و ۵۰۰ ماده) در 
سل ۵۱ گسترش یافت. این اندازه جمعیت برای ٩‏ نسل بعدی (تا 
ذزسل ۶۰) ثابت ماند. تنها افراد سل ۵۱ دارای رکورد فنوتیپی بودند و 
گروه مرجع (گروهی که از اطلاعات آنها جهت برآورد اثرات ذشانگری 
استفاده می شود) را تشکیل می‌دادند و افراد سایر نسل‌ها (۵۲ تا ۶۰) 
فا قد اطلاعات فنوتیپی بود ند و گروه تأّیید (گروهی که فاقد رکورد 
فنوتیپی بوده و ارزش اصلاحی زنومی آنها تنها بر اساس اطلاعات 
ژنوتیپی برآورد می‌شسود) را به وجود آوردند. سبیه‌سازی داده‌ها با 
استفاده از نرم‌افزار حطاعحصن (۱۵) انجام شد. 

جهت سنجش صحت شبیه‌سازی از روش بررسی روند عدم تعادل 
پیوستگی در هر نسل استفاده شد. برای این منظور مقدار مورد انتظار 
عدم تعادل پیوستگی در هر نسل بر اساس فرمول اسوید (۱۷) محاسبه 
شد: رابطه ۱ 
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در رابطه فوق (۳)۳7 میزان مورد انتظار عدم تعادل پیوستگی در 
نسل 0ام» مل] اندازه موّثر جمعیت و > میزان نرخ نوترکیبی بر | ساس 
مورگان هستند. 
شبیه سازی ژنوم 

در این مطالعه ژنومی شامل سه کروموزوم هر کدام به طول یک 
مورگان در نظر گرفته شد. روی هر کروموزوم سه سطح تراکم 
چندشکلی تک نوکلئوتیدی دو آللی شامل (۸۵۰۰ ۱۰۰۰ و ۱۵۰۰) در 
مورگان با مقادیر مختلف 941 شامل ۵۰ (پایین) ۱۰۰ (متوسط) و ۱۵۰ 
(بالا) 40 با فوااصل یکسان نشانگری #سترده شدند. همچنین جهت 
مطالعه تأثیر ورائت‌پذیری و نحوه توزیع واریانس 9/1 بر صحت برآورد 
ارزش‌های اصلاحی ژنومی از سه صفت گوسفند با ورائت‌پذیری پایین 
(۰/۱)» متو سط (۰/۳) و بالا (۰/۵) و دو نوع توزیع واریانس 901 شامل 
توزیع یکنواخت و گاما (8-۰/۴ و 0-۱/۶۶) (۷) در شبیه‌سازی استفاده 
شدند. هم نشانگرها و هم ها دو آللی بودند. 

این اسان ۱۸ شتارتوی کر کییی یف کل ای ۳ ود 
برای تعداد 411» ۲ سطح برای نحوه توزیع واریانس 411 و ۳ سطح 
برای تراکم ذ شانگرهای 900 می‌با شد. همچنین سطوح ورائت‌پذیری 
در هر کدام از سناریوها تغییر داده شد تا تأثیر آنها بر صحت برآورد 
ارزش‌های اصلاحی ژنومی بررسی شود. 


پیش‌بینی ارزش های اصلاحی ژنومی 
روش 681/۳ 

روش ,]۳ به کار رفته در این مطالعه بر اساس یک مدل 
مختلط ساده طراحی شد و فرض گردید که در آن کلیه چندشکلی‌های 
تک نوکلئوتیدی بر صفت مورد مطالعه موّثر می‌باشند. برای به کار 
بردن این روش از نرم‌افزار اجع۸۹1 (۶) و مدل زیر استفاده شد: 


مقایسه روش‌های 68111۴ و بیزی در برآورد ارزش‌های اصلاحی ژنومی با معماری‌های... ۰ ۲۳۳ 


رابطه ۲ 21٩‏ + برللا < 7« 


جدول ۱- ساختار جمعیتی و متغیرهای استفاده شده در شبیه‌سازی 


وه امه تهج ممتاملنته 4ص منک ممت)‌هلنام۳0 -1 ماماه 1 
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۰ رأس دام (۵۰ رآس نر و ۵۰ رأس ماده) 


۰ رأس دام (۵۰۰ رأس نر و ۵۰۰ رأس ماده) 
تابع مکان‌یابی هالدان 


60 


10 


که در آن لا بردار ارزش های فنوتیپی» ۸ میانگین فنوتیپی» 9 
تعداد رکورد هاء مب بردار یکا به ابعاد 1<*0 2 بردار در برگیر نده 
ارزش های اصلاحی ژنومی. 7 ماتریس ضرا یب ارتباط دهنده 
مشساهدات به بردار ارزش‌های اصلاحی و ه بردار اثرات تصادفی 
باقیمانده می باشد. اثرات ژنتیک افزایشی چندشکلی‌های تک 
نوکلتوتیدی (ع) دارای توزیع نرمال (م1۷)0,6 فرض شد که 27 
ماتریس روابط خویشاوندی برای همه چندشکلی‌های تک نوکلئوتیدی 
می‌باشد که و6 انحراف معیار اثرات ژنتیکی ها می‌باشد. باید 
توجه دا شت که 2 بر | ساس احتمالات 135 بین هر دو دام محا سبه 


۱: 

اندازه موّثر جمعیت 

(۳۲9) 5126 طمتاه000 ۳۲۲۵۵۱۷۵ 
طول ژنوم 

طاعمع1 معمصعت) 

تعداد چندشکلی‌های تک نوکائوتیدی 
۱۱۱ 
تعداد ب] 001 

1 ها 

توزیع اثر با ]0 

عصمتا ماوت واهعه اه 

ورائت پذیری صفت 

(2 7 

جمعیت مرجع 

۶ عصلطل۲2 1[ 

۲۷۵11021102 6 

تعداد نسل در جمعیت عدم تعادل پیوستگی 
0ص م1 فممتاهمصوهع ۵۶ بعماصی 
هام۵0 مامتان وع وت 

ندازه جمعیت عدم تعادل پیوستگی در هر نسل 
۱۹ 

ممتاهتعصعم هه میاتتماتانن وه ول 
اندازه جمعیت 

۳00۱9124100 6۵ 

نوترکیبی 

1۴۵6001210۲ 

تعداد نسل جمعیت 


0 0 000۱۱2 0۶ تمطاصت 


۳ 


فصمتاتاه م۲6 0۶ مایا 


شده است و برای تمام دام‌های دارای ژنوتیپ (دام‌های جمعیت مرجع 
و جمعیت تأیید) تعمیم یافته است. برای محاسبه مجموع روابط آللی 
برای پیش‌بینی اثرات افزایشی چندشکلی‌های تک نوکلئوتیدی از 
مدل‌های مختلط زیر استفاده شد: 


رابطه ۳ 


۳ ِ 7 12 1 1 
7 1۵+ 72 ب71 


۲ _ ۳ 2 ۹ 
که وان یی 0 و ارمانزیشن همانی ام وا 
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روش بیز ۸ 
در همه روش‌های بیزی از مدل آماری زیر استفاده شد: 
رابطه ۴ 6 + :9ز ی + ملم  <‏ 


که در آن ز بردار ارزش های فنوتیپی» ۸ میانگین فنوتیپی» 9 
تعداد رکوردهاء بر بردار یکا به ابعاد 1<0, :6 ماتریس ژنوتیپ‌ها 
برای 1 امین چند شکلی تک نوکلئوتیدی, رو اثر جایگزینی اللی [ امین 
چندشکلی تک نو کلئوتیدی و 6 بردار اثرات تصادفی باقیمانده می‌باشد. 

در روش ۳2765۸ همانند روش نانآ فرض می‌شود که 
تمامی ٩۳‏ ها دارای اثره حتی جزئی بوده و تعدادی از ٩۳‏ ها در 
دارند. اثرات ها از توزیع نرمال با واریانس جداگانه برای هر ٩۳/۳‏ 
از توزیع کای دو معکوس» نمونه‌گیری می‌شود (۱۰). برای به کار بردن 
این روش از معادلات زیر استفاده شد: 

2 2 2 2 
۵ رورم ۵ (تو تاره ۶ 
21 
و2 (1-7) 2 
(2 - )م65 و 
سس و 
۳ ۳ ۷" 
که در آن" دلالت بر درجه آزادی و 5 مربوط به پارامتر معیار است. 


روش بیز ظ 

در اين روش اکثر ٩(8‏ های موجود در منطقه ژنومی 011ای 
نداشته و بنابراین اثرات آنها برابر با صفر می‌باشده در حالی‌که تعداد 
اندکی از آنها (-۱) در عدم تعادل پیوستگی با ب011 قرار داسته و 
دارای اثر می‌با شند. در نتیجه اثرات غیر صفر در عدم تعادل پیو ستگی 
بالایی با 011 قرار دارند. بیز ط تحت معادلات زیر اجرا شد: 


رابطه ۶ 

0-۵ نامهتم (قور) 2 َ( مر 
۲ ۲۵۵۵۵ ۱۳ 0 نّ 

روش بیز ) 


فزاروش بر 6 هن نی نی سیگ که توت میتی از 
چند شکلی‌های تک نوکلئوتیدی (7) بر صفت اثر دارند. در واقع بیز » 
شکلی از بیز 13 است. فقط به جای استفاده از توزیع پیشین ] برای 
تأثیرات نشانگری, از توزیع نرمال استفاده می‌کند. برای به کار بردن 
اين روش از مدل زیر استفاده گردید: رابطه ۷ 

۱ 6 ۲ :2:9 + م1 < ۲ 

که در آن ۷ بردار ارزش‌های فنوتیپی» ۸ میانگین فنوتیپی» 1 تعداد 
رکوردهاء بر بردار یکا به ابعاد ۰1*0 .2 بردار ژنوتیپ ها برای امین 
چندشکلی تک نوکلوتیدی, ,0 اثر جایگزینی آللی‌امین چندشکلی 


تک نو کلئوتیدی و » بردار اثرات باقیمانده‌است که دارای توزیع نرمال 


مس اف 


روش بیز 1۲۸0 

در روش بیز 12850 از تابع نمائی دو گانه برای توزیع اثرات :0011 
و ٩(/۳‏ تحت مدل بیزی استفاده شد: 
رابطه ۸ (۷)۵,۱0,۰2۲۶ ] < (جّه,۸|رط) 5 

وه 
ز6 7 1 ۳ 

که در معادله فوق .2 ناشناخته فرض شده و به آن توزیع گاما با 
پارامتر شکل .2 و مقیاس ۲ اختصاص داده می‌شود. 

برای محا سبه ارزش | صلاحی واقعی حیوانات از رابطه ٩‏ | ستفاده 
شد: 

رابطه ۹ رفازن و-زرظ كت 1۳9۷ 

که در معادله فوق 7817 ارزش اصلاحی واقعی افراد» 0 تعداد 
بآ های موثر بر صفت. رز تعداد آلل‌های ۱ که فرد ز در لوکوس [ٌ 
حمل می کند و 2 اثر [امین 01 می‌باشد. در این مطالعه همه 
۰ رأس دام نسل ۵۱ به عنوان جمعیت مرجع به کار گرفته شدند و 
از فنوتیپ دام‌ها در اين ذ سل برای برآورد ارزش‌های | صلاحی ژنومی 
در جمعیت تأیید استفاده شسد. ارزش اصلاحی ژنومی برآورد شسده 
ارزش اصلاحی ژنومی افراد گروه تأید) از مجموع آثار نشسانگرها با 
توجه به ژنونیپ آنها به صورت زیر محاسبه شد: 

اب 8 2 658۷ 

که در آن 26 ماتریس طرح می‌باشد که به ژنوتیپ دام‌ها اشاره 
دارد و 2 بردار اثرات برآورد شده برای چندشکلی‌های تک نوکلئوتیدی 
می‌باشد. 
ارزش‌های | صلاحی واقعی و ارزش‌های | صلاحی ژنومی برآورد شده 
(۱0۷,۵00۷) محاسبه شد. هر آزمایش ۱۰ بار تکرار و سپس میانگین 

در اجرای روش های بیزی از روش نمو نه‌برداری گیبس و نرم 
افزار > ورژن ۲-۱۰۱ (۱)۱۳ستفاده شد. برای هر تحلیل یک زنجیره 
مارکوف مونت کارلو تا رسیدن به همگرایی به طول ۲۱۰ هزار چرخه 
شامل ۱۰ هزار چرخه اولیه جهت گرم شدن و ۲۰۰ هزار چرخه اصلی 
به کار رفت که در هر ۵ دور یک بار نتایج ذخیره شد و در پایان نتایج. 
۰ هزار چرخه ذخیره گردید. 


میانگین عدم تعادل پیو ستگی (*) افراد سل آخر جمعیت (ذسل 
۰ م) ۰/۰۱۱ ۰/۱۹۱ بدست آمد. این در حالی است که مقدار مورد 
انتظار آن توسط فرمول اسوید (۱۶) ۰/۲۱۹ برآورد گردید. این مقدار 
نشان می‌دهد که به طور میانگین بیش از ۸۷ درصد از حداکثر مورد 


مقایسه روش‌های 68111۴ و بیزی در برآورد ارزش‌های اصلاحی ژنومی با معماری‌های... ۰ ۲۴۳۵ 


انتظار محقق گردیده است که میزان مطلوبی است. میانگین عدم 
تعادل پیوستگی با مطالعات ویلیامسن و همکاران (۱۷) مطابقت دارد. 
تأثیر تعداد 011 بر صحت برآورد ارزش‌های اصلاحی ژنومی 
صحت ارزش‌های اصلاحی ژنومی برآورد شده با سه سطح تعداد 
01 و ورات‌پذیری در جدول ۲ ارائه شده است. در هر پنج روش 
آماری با افزایش تعداد 91ها از ۵۰ به ۱۵۰ صحت ارزش‌های ژنومی 


بالاترین صحت برآورد مربوط به روش بیز 9 و کمترین مقدار توسط 


روش 801110 حاصل شد. این نتایج با نتایج دیتویلر و همکاران (۵) که 
گزارش کرده‌اند با افزايش تعداد 41 موثر بر صفت» صحت برآوردهای 
ارزش‌های اصلاحی ژنومی کاهش می‌یابد. مطابقت دارد. 

دلیل اين امر را می‌توان به توزیع مقدار محدودی واریانس ژنتیکی 
بر تعداد زیادی 41 دانست که در نتیجه سهم هر 0و در ارزش ژنتیکی 
کل کاهش يافته است و قدرت مدل‌ها در برآورد تأثیرات کم خواهد 


جدول ۲- صحت ارزش‌های اصلاحی ژنومی برآورد شده با سه سطح تعداد -011) و ورائت‌پذیری 
تدم اک مصصعء صا راتاتطاماتعط فصح مته‌طاصصه با 00 معط جوا (عو) بچمهساهعه متعطمجمع 12060)وه 1۳6 -2 12016 
مدل آماری 
90061 21501021 


۱ 1-0 ۳2۷7۶ 0 
0.757)0.01( ۰ 07560.02( ۰۰. 075760.02( 
0743)0.01( ۰۰. 0741)0:02( ۰... 07400.02( 
0.658)0.01( 0.652)00.02( 0.65200.02( 
0.776)0.01( ۰ 077000:02 ۰.۰. 076760.02( 
0768)0.01( ۰۰. 076400:02 ۰۰. 07620.02( 
0:657)0.01( ۰۰... 0:652)0:02( ۰۰. 06500.02( 
0.840)0.01( ۰۰ 0:838)0.01( ۰ 08360010 
082100:01 ۰. 0818001( ۰. 0817000 
0743)0.01( ۰۰. 0739001 ۰... 07350010 


مطالعات در این زمینه نشسان داده است که وقتی تعداد 011های 
موّثر بر صفت زیاد باشد روش 810 مشابه و يا حتی بهتر از 
روش‌های بیزی عمل می‌کند (۵ و ۱5). با وجود این» به دلیل اينکه 
روش‌های بیزی معرفی بهتری از نحوه توزیع اثرات 011 بر صفات» 
نسبت به روش 03۳ به صورت پیش‌فرض ارائه می‌دهد در تعداد 
اندک 0 برتر عمل می‌کنند. 

همچنین صحت برآورد ارزش‌های اصلاحی ژنومی صفات با 
ورائت‌پذیری بالا ذ سبت به ورائت‌پذیری پایین بي شتر ا ست. دلیل این 
مسئله واضح است. زیرا هرچه ورائت‌پذیری صفت بیشتر باشد, فنوتیپ 
فرد به ارزش ژنتیکی فرد نزدیک‌تر بوده و در نتیجه آثار نشانگرها و به 
دنبال آن ارزش اصلاحی ژنومی افراد به طور صحیح‌تر برآورد می‌شود 


اد ۰ وش 

1-۸ طلاظات ‏ ولا وتانطمانه۲ 
(0.685)0.01 (0.67900.01 50 

01 100 0.673)0.01( 0.681)0.01( 
150 0.66900.01( 0.672)0.01( 
50 0.689)0.01( 0.695)0.01( 

03 100 0.681)0.01( 0.687)0.01( 
150 0.678)00.01( 0.67500.01( 
50 0.756)0.01( 0.75900.01( 

05 100 0.7270.01( 0.73400.01( 
150 0.715)0.01( 0.727)0.01( 

.)۱۸( 


تأثیر تراکم نشانگر بر صحت پیش‌بینی ارزش‌های اصلاحی 
ژنومی 

صحت ارزش‌های اصلاحی ژنومی برآورد شده با سه سطح تعداد 
بل مراکم تشانگری در و۲۰ راک شده ایبه: .با افرآن تراک 
ذشانگرها صحت پیش‌بینی تو سط هر پنج روش آماری افزايش یافت. 
نتایج حاصل از این پژوهش با نتایج تحقیقات کالوس و همکاران (۲) 
و مویر (۱۱) مطابقت دا شت. بیشترین صحت ارزش | صلاحی ژنومی 
توسط روش بیز 0 در تراکم نشسانگری ۱۵۰۰ و تعداد ۵۰ ,]01 به 


دست آمد. 


۳۴۶ نشریه پژوهشهای علوم دامی ایران جلد ۰۱۲ شماره ۲ تابستان ۱۳۹۹ 


جدول ۳- صحت ارزش‌های اصلاحی ژنومی برآورد شده با سه سطح تعداد 061 و تراکم نشانگری 
٩۳۱۳ 1005107 6‏ ره فتعحاصصباه ب 01 عمط ت۲۵ (مو) پرمهتنهمج متمصمصمع 0عامحصطتاوه 3-1 م1201 


مدل آماری 
9061 1و1و2 
تراکم نشانگری ‏ تعداد 011 
0ص ۳2۲۵5 1-0 و۳26 ۳2۸ ات «اتعصعل ٩۳۲۲‏ ۱۹/۹ 
 0:68700.01( ۰ 0:76900:01( ۰ 0762)0.02( 0.760)0.02(‏ (0:685)0.01 500 
 0.77200.01( 0.768)00.02( 0.763)0.02(‏ (0.69300.01 ۰ (0.68800.01 1000 50 
 0.78800.01( 0.780)00.02( 0.779)0.02(‏ (0.70500.01 (0.694)0.01 1500 
(0.756)0.02 (0.758)0.02 (0.76100.01 ۰ (0.67700.01 (0.67600.01 500 
(0.759)0.02 (0.761)0.02 (0.76500.01 ۰ (0.68600.01 ۰ (0.67900.01 1000 100 
(0.771)0.02 (0:77300.02 ۰ (0:7780.01 ۰ (0:69900.01 ۰ (0.687)0.01 1500 
(0.641)0.01 (0:648)0.01 ۰ (0:64900:01 ۰ (0:66000.01 ۰ (0:669)0.01 500 
(0/648)0.01 06510010 .۰ 065500:01 ۰ (0:67300:01 ۰ (0:675)0:01 1000 130 
(0.661)0.01 06656001 .۰ (0:6680:01 .۰ (0:68900:01 .۰ (0:692)0:01 1500 


سولبرگ و همکاران (۱۶) گزارش کردند که با افزايش ذ شانگرها 
از ۱۰۰ به ۸۰۰ نشانگر در هر مورگان صحت ارزیابی ژنومی از ۶۹ 
درصد به ۸۶ درصد افزایش یافت. افزايش تعداد نشسانگرها باعث 
فزایش میزان عدم تعادل پیو ستگی بین ژن‌ها و شانگرها می شود و 
بدین ترتیب صحت ارزیابی ژنومی را افزایش می‌دهد. 


تأثیر ورائت‌پذیری بر صحت پیش‌بینی ارزش‌های اصلاحی 


ورائت‌پذیری و تراکم ننشانگری در جدول ۴ ارائه سده است. نتایج 
نشان داد که افزايش تراکم نشانگری در سطوح ورائت‌پذیری پایین 
(۰/۱) و بالا (۰/۵)» سبب افزايش صحت ارزش اصلاحی ژنومی 
می‌شسود اما در صفات با ورائت‌پذیری متوسط (۰/۲) افزایش تعداد 
نشانگرها بر صحت ارزیابی ژنومی تأثیری نداشت. بطوری که بیشترین 
صحت پیش‌بینی ارزش اصلاحی ژنومی در ورات‌پذیری ۰/۵ با تراکم 
ذشانگری ۱۵۰۰ تو سط روش بیز 3 و کمترین صحت. تو سط روش 


و لاه در ورات‌پذیری ۰/۱ با تراکم ذشانگری ۵۰۰ به د ست آمد. 
ر ۱ 
۱ شب اين نتایج با پژوهش‌های پیشین مطابقت داشت (۱۴ و ۷ 
یط ین نتایج با پژوهش‌های پیشین (۱۴ و ۷) 
جدول - صحت ارزش‌های اصلاحی ژنومی برآورد شده با سه سطح تراکم نشانگری و ورائت‌پذیری 
هنعط مه راتعصعل طا5 عععطا 1۵۲ (0و) پرمه‌تنهم2 متصمصعم 60امحصتاوع 12۳0164-16 


مدل آماری 
۱۱ 
تراکم نشانگری. ورائت پذیری 
۵0 ۳2۷۵9 اقت ۱ ۱ ۱ ناه حاتعدعل ٩(۱۲‏ «انازطامات۲۵ 
(0.73900.02 (0.740)0.02 ۰۰ (0:74200.01 .۰۰ (0:68000.01 .۰ (0:67300.01 50 
(0.74400.02 0751002 . 0759000 0699000۵ .۰ (0:6820:01 1000 01 
(0.754)00.02 (0758)0.02 ... 07630010 070700۵ 069500:01 1500 
(0.76000.02 (0.760)0.02 ۰۰ (0:76400.01 .۰ (0:69500.01 .۰۰ (0:692)0.01 50 
(0.76000.02 0760002 ...07650001 069500 ...۰ (0:69300:01 1000 03 
(0.761)00.02 0761002 0766000 069600۵ ...۰ (0:693)0:01 1500 
(0.77100.01 (0771)0.01 ۰۰ 077900:01 .۰۰ (0:71900.01 .۰ (0.715)0.01 50 
0.77900.01 0781000 ...۰ 07860010 .۰ 07290010 .۰۰ (0:72200:01 1000 05 
0.785)00.01 (0788)0.01 .. 07910010 074200:01 ...۰ (0:73800:01 1500 


موویسن و همکاران (۱۰) و گودارد (۷) یکی از عوامل اصلی موثر 
بر صحت انتخاب ژنومی را ورات‌پذیری صفت می‌دانند. بي شتر بودن 


فنوتیبی, به معنای بیشتر بودن نقش ژن‌های با بیان افزایشی در ایجاد 
پراکنش در صفت است که باعث برآورد صحیح‌تر تأثیرات ذ شانگرها 


ورائت‌پذیری» یعنی بيٍ شتر بودن ذ سبت واریانس افزاي شی به واریانس می شوند. صحت برآورد ارززش اصلاحی ژنومی صفات با ورائت‌پذیری 


مقایسه روش‌های 68111۴ و بیزی در برآورد ارزش‌های اصلاحی ژنومی با معماری‌های... ۰ ۲۳۷ 


بالا نسبت به ورائت‌پذیری پایین بیشتر است. دلیل اين مسئله واضح 
است زیرا هرچه ورائت‌پذیری صفت بیشتر باشد فنوتیپ فرد به ارزش 
ژنتیکی فرد نزدیک‌تر بوده و در نتیحه آثار نشانگرها و به دنبال آن 
ارزش اصلاحی ژنومی افراد به طور صحیح‌تر برآورد می‌شود (۱۸ و ۷). 


نتایج تأثیر واربانس :011 بر صحت برآورد ارزش‌های 
اصلاحی ژنومی 

صحت ارزش‌های اصلاحی ژنومی در تراکم نشانگری و تعداد 
17 یکسان با توزیع واریانس یکنواخت و گاما در جدول ۵ ارائه شده 
است. صحت ارزش‌های اصلاحی در حالت توزیع گاما تأثیر ژنی 
بهتری ذ سبت به توزیع یکنواخت دا شت. همچنین. هنگامی که توزیع 
تأثیر ژنی گاما است روش‌های بیزی در مقایسه با روش دااداع از 
صحت ارزش ژنومی بالاتری برخوردار هستند. در بین روش‌های 
بیزی» روش بیزط بیز) و بیز۸۹80] در حالت توزیع گاما تأثیرات 
ژنی بهترین عملکرد را نان داد ند اما تفاوت آذها از نظر آماری 
معنی‌دار نبود (0۰/۰۵). نتایج حاصل از این تحقیق با گزارش‌های 


موویسن و همکاران (۱۰)» گودارد (۷)» سولبرگ و همکاران (۱۶) و 
دیتویلر و همکاران (۵) مطابقت داشت. 

دلیل اين امر دو عامل می‌تواند باشسد: نخست آنکه توزیع پیش 
فرض روش های بیزی با توزیع گامای وار یانس ژن های عمده اثر 
هماهنگی دارند که باعث برآورد صحیحتر اثرات چندشکلی‌های تک 
نوکلئوتیدی می‌گردد. همانطور که گودارد (۷) گزارش کرده است توزیع 
پیش فرض گاما در روش های بیزی برآوردهای صحیح‌تری ایجاد 
می‌کند که با نتایج این مطالعه هم‌خوانی دارد. دلیل دوم این است که 
توزیع گاما دارای تعدادی ژن‌های عمده اثر با اثرات بزرگ‌تر ذ سبت به 
توزیع یکنواخت است. در نتیجه برآوردهای ارزش‌های اصلاحی ژنومی 
در صفاتی با تعدادی ژن‌های عمده اثر با اثرات بزرگتر صحیح‌تر است. 

یکی از ضعف‌های روش 201 این است که برای همه نشانگرها 
سهم یکسانی از واریانس را در پیش‌بینی ارزش‌های | صلاحی در نظر 
می‌گیرد. این در حالی است که در روش‌های بیزی بر حسب توزیع 
پیشین» وزن‌های متفاوتی به نشانگرها اختصاص داده می‌شود. 


جدول ۵- صحت ارزش‌های اصلاحی نومی در تراکم نشانگری و تعداد له یکسان با توزیع وریانس یکنواخت و کام 
مصمتنه۷ مصصصصعع 0صه متمگتصا طز۳ تاتعجعل ٩۲۲‏ 4صع وتمجامصنه ]00۲ محصصعع صا تدممتتامه متحعمصعع 0عتمصصتاوع م19 -5 12016 


ممتاباطاتتاونل 
مدل آماری 
0061 ۹2115)1621 
ورائت پذیری توزیع واربانس 
0۷2 
0 ۳2۷ ۱۵ ۱ ۱2 نان هط 
(0:712)0.02 (0:715)0.02 ۰۰ 071900:01 .۰ (0:677)0:01 .۰ (0:67400.01 01 ۳1 
(0.726)0.02 0730002 ...073400 0685001 .۰ (0:68300:01 03 ِ 
(0.796)0.02 080100۵ ...۰ 080200 ...۰ 0740001 ...۰ (0733)00:01 05 سولتسنا 
(0.753)00.02 (0.758)00.02 (0.760)00.01 (0.691)00.01 (0.684)00.01 01 اما 
(0.776)0.02 077600۵ 0779000 070000 (0:695)0:01 03 
(0.841)0.02 0842002 ...0845000 0749001 ...۰ (073900:01 05 شسست 
نتیجه گیری کلی مطلوب‌تری ارائه دهد. لیکن چنانچه صفت تحت تأثیر تعداد زیادی 


نتایج این تحقیق نشان داد که عواملی نظیر تراکم نشانگرهاء تعداد 
1 توزیع واریانس 011 و ورات‌پذیری صفت بر میزان صحت 
برآورد ارزش های اصلاحی ژنومی موّثر هستند. در صفات با 
ورائت‌پذیری بالا هر چه تراکم نشانگرها بیشستر و تعداد 0 کمتره 
مطالعات ژنومی چنانچه صفتی تحت تأثیر تعداد کمی 0و قرار دا شته 
باشد. برآورد ارزش‌های اصلاحی توسط روش بیزه می‌تواند نتیجه 


ست. | ستفاده از هر پنج روش آماری «داااع. بیزه بیزط بیزن و 
بیز 12550 مزیت یکسانی خواهد داشت. حتی در بعضی موارد روش 
«دااع عملکرد بهتری نس بت به چ هار روش دیگر دارد. تراکم 
ذشانگرها بر صحت ارزیابی ژنومی در صفات با ورائت پذیری متو سط 
تأثیری نداشت و از آنجا که بیشتر صفات اقتصادی در نژادهای بومی 
گو سفند از د سته صفات با ورائت‌پذیری متو سط ه ستند. می‌توان از 
۰ تا ۱۰۰۰ نش‌انگر برای برآورد ارزش‌های اصلاحی ژنومی در 
برنامه‌های شبیه‌سازی استفاده نمود. 
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۷ 20101016 130102500 ر۱۵۷۵[5 01110ه)۵11 (0.5) طونط 4ج (0.1) ۱0 26 پچاتعصعل تفج مصتفههتمصا اقط 530۱7۵۵ 
۶ مدع ,صمامنلدنه عتعصمصعع ۵۶ مدمه مط امعله امظ 010 فاتهتا (0.3) بوتلتماه)اتعط ماقعلمصه 2 مطلوهعمصا غیاط 
۰ متا ما قامعره مصمع باه وم۳۵۲10 ممتباطاتتادنل مصصصصوع مط مر فمیتاه۷ عصتلعع۳ متمتمصوع 


طعمطمهط۱۱۵ ۵۶ ازوهب نصا روهم‌تنمفم؟ امتیتدا مه مسالنه‌ته۸ ۵۴ ولنهد۲ رقع‌مممهه تمصتمظ ۵۶ ,12006 ریگ۳0ظ 1-۸9500126 
۰ ,۸۲02011 ,۸۲020111 

20 رملتطاه1 1۵02 ۵۶ نویه نصا روه‌مممزمو امسانه‌نهد ۶ه وانمدظ روم‌مممتمه امصصتمظ ۵۶ 1۳۵۵۲۰ رله‌اهن0متم 2-۳ 
(صمه.ممطهبم6) زا تمطو0عبمو ۲۵2۵ :1نممطر! عمطانته ممتلم0موعترمن) -۲) 

101: 10.22067/1[257- 1.2 


۱۲۳۵۰ نشریه پژوهشهای علوم دامی ایران جلد ۰۱۲ شماره ۲ تابستان ۱۳۹۹ 


0ص )0ع1ه ب1 0۵ ط0ز)ناحا زاون روته‌جاحصناط با1 ی ووتتعصهمل تمانهوه قح طمیاع قمامه) اقط ۵0«مدو فالناوعم م1 صمتقوت مومت 
,کل 1110حا2)تبعط طونط ما .ما۷ ممتلععته متمصمصعع ۵۶ بممتنهعه معط مصتامصتاوه صا وتا0عه متع۷ باتازماهاتبمط نم 
هصتل0۳66 متعمصمع عصتامصتایی آه بمعسهعه 0و2ع07ص1 وه ص162 روته‌طامصتاه 021 م1 قمع بوتعصمه ومآنعمط تمطم‌نط فطل 
7 عمباله۷ عصتلمعه ۵۶ ممتاقتاده رصا1 ۵ ۵۶ تمجاحصباه ألقصصری 2 بوطا 0متامع21 و1 نا مط 1 رومتنتاد متصطمصمع ما .عمناله۷ 
متصمصهع ۵۶ پمدتساهمه مطا اعد امط 010 فمتاتفصع بملیها ,البافعر مامام1270 وعمحظ ج ۱۷۸۱6۱4 صو 0طاعحظ ظ فرظ 
وممطو گم کم50601 21170 ط1 صاتهتا متممجممع فط ۵۶ )ومصظ ومصی 24 وزاتااماه‌اتهط مقهمصهط گم فاتقط ما ممتامتآهه 
0۳۵22 مرمنج آباصزی ما فمبله۷ ممتلهعره ماقحصتاوع م6 0عون1 عم وم وتماتهجه 1000 مه 500 ات201 تعط 2۵067246 


۰ تاباظ تمصع رمففضصا مهتفعروظ ر2۵9ظ ,فورظ رففع27ظ :15۵0۱۷۵۲۵۵ 


